Spectrum of the Vortex Bound States of the Dirac and Schrodinger Hamiltonian in 

the presence of Superconducting Gaps 
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We investigate the vortex bound states both Schrodinger and Dirac Hamiltonian with the s-wave 
superconducting pairing gap by solving the mean-field Bogoliubov-de-Gennes equations. The exact 
vortex bound states spectrum is numerically determined by the integration method, and also ac- 
companied by the quasi-classical analysis. It is found that the bound state energies is proportional 
to the vortex angular momentum when the chemical potential is large enough. By applying the 
external magnetic field, the vortex bound state energies of the Dirac Hamiltonian are almost un- 
changed; whereas the energy shift of the Schrodinger Hamiltonian is proportional to the magnetic 
field. These qualitative differences may serve as an indirect evidence of the existence of Majorana 
fermions in which the zero mode exists in the case of the Dirac Hamiltonian only. 

PACS numbers: 71.27.+a, 74.45,+c, 71.10.Pm 



I. INTRODUCTION 

The zero mode (Majorana fermion excitation) attracts 
lots of investigation due to its non-Abelian statistics 0- 
0] and possible application to the fault tolerant quantum 
processing [3, Possible candidates to support Majo- 
rana fermions are p + ip superconductors @, S 0] > P + W 
supcrfluids in cold atoms |8l4ldl|. With the proximity ef- 
fect between the s-wave supe rcondu ctor and the strong 
topological insulator surface 
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chiral Majorana 

fermions could be created as edge states [HI]. Besides 
that, Majorana fermions could also be realized in semi- 
conductor with spin-orbit coupling 16l- 18|. 

To verify the existence of Majorana fermion, Law et al. 
proposed the tunneling experiments to probe the chiral 
Majorana fermion at the interface between a supercon- 
ductor and the surface of a topological insulator, in which 
the Maj orana fermions induce resonant Andreev reflec- 
tion 0. 

In this paper, we investigate and compare the vortex 
bound states of both the Dirac and Schrodinger Hamilto- 
nian with the s-wave superconducting gap by solving the 
mean-field Bogoliubov-de-Gennes equations. The bound 
states are numerically solved by the integration method 
as well as the quasi-classical analysis. Their differences 
could be served as an indirect verification of the existence 
of the Majorana fermion at the surface of the topologi- 
cal insulator with an induced superconducting gap due 
to the proximity effect. 



II. FORMULATION OF VORTEX BOUND 
STATES 



Our formulation is based on the mean-field 
Bogoliubov-de-Gennes (BdG) equations for the quasi- 
particle. In terms of Nambu indices, the mean-field 



Hamiltonian becomes 
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where = (c t ,c ; ) T and $ = (c\,c\) T . H { Q D ' S) 
paper could be the Dirac Hamiltonian, Hq D ^ = v-pp- a — 
/i — ha z , and the Schrodinger Hamiltonian, Hq = £— — 
fi— ha z . v-p denotes the Fermi velocity, \i is the chemical 
potential, h represents the Zeeman coupling due to the 
external magnetic field acting along the z-direction, and 
A(r) is the superconducting gap function. 

The quasi-particle states of the above Hamiltonian in 
Eq.([T]) are given by 
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Considering an isolated vortex carrying one flux quan- 
tum, i.e. A(f) = A(r)e 1 ^, with the eigenvectors of the 
form, u n (r) = (u^ n {r)e m * , H . n {r)e^ n+1 ^) and v„(f) = 
(w ti „(r)e™*,'y i:n (r)e i (™- 1 ^). The eigenproblem for the 
Dirac Hamiltonian in radial coordinate becomes 
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where 

-- 

and u n (r) = (u ti „(r), u^ n (r)Y , v n (r) 

{vt, n (r),vi, n (r)) T . 

For the Schrodinger Hamiltonian, the eigenproblem 
can be further block-diagonalized into 

u tA r ) \ = E ( "t.nO) 

(5) 
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Eqs.© and ([6]) are not independent to each other. In 
fact, their complete sets of eigenvectors are the same. 

The size of the vortex core is characterized by the co- 
herent length £. For the Dirac Hamiltonian, it should be 
of the order of %vp / Aq , the unique length scale of the sys- 
tem, where Ao is the superconducting gap far away from 
the vortex core. However, for the Schrodinger Hamilto- 
nian, there are two length scales kp 1 and Hvf/ Ao. £ is of 
the order of Hvf/A when /i/A ~> 1. As /i/Aq becomes 
smaller, £ is going to shrink to another length scale, fcp 1 
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26]. When fJ,/A turns to be negative, a new length 
scale due to the bosonic molecules emerges. 

We adopt £ = fwy / Aq for the Dirac Hamiltonian, and 
£ = he/ Aq for the Schrodinger Hamiltonian in which c 
is a constant independent of the chemical potential [i. 
Further we assume the form of the gap amplitude 



A(r) = A tanh(-) 



(8) 



The BdG equations, Eqs.© and ([5]), are solved by 
the integration method. Imposing that the wavefunction 
aymptotically approaches to zero far away from the vor- 
tex core, the eigenenergies can be obtained by integrating 
the differential equation from the infinity to the vortex 
center. The advantage of the integration method is that 
we can solve for the numerically exact (that is, system 
size in the thermodynamical limit) bound state energies 
once the form of the gap amplitude as in Eq. jSJ is given. 

The chemical potential fi in the BdG equations for 
the Schrodinger Hamiltonian can be positive or nega- 
tive, arbitrary small or large. The condition of the 
quasi-classical approximation, fcp£ 3> 1 (equivalent to 
/j,/Aq ^> 1), may not be valid in the whole range of the 
chemical potential. Eq.© is solved directly by the inte- 
gration method. We do not assume the Andreev (quasi- 
classical) approximation 23- 25j|. 



III. NUMERICAL RESULTS 

We first consider the case without external magnetic 
field. The vortex bound state energies of the Dirac 
Hamiltonian are solved for different chemical potential. 



The vortex bound state energies of the Dirac Hamil- 
tonian as function of chemical potential /i is shown in 
FigfT] We do not show the negative n because E_ n — 
—E n which is the consequence of time- reversal symmetry. 
That is, E — > —E when n — > —n (reverse the angular mo- 
mentum). Besides that, E — > —E when \x — > —fi, which 
is the particle-hole symmetry particularly for the Dirac 
Hamiltonian. 

For n = 0, the zero mode corresponds to the Majorana 
fermions. For non-zero angular momentum, n > 0, 3> 
Ao, we found the mid-gap state near zero energy. As /i 
approaches to zero, the mid-gap state gradually merges 
to the continuum of energy E — Aq. Especially, near zero 
chemical potential, fj, ~ 0, all the vortex bound states 
except the Majorana fermions state, are located near the 
edge. 

We also plot the bound state energy as a function of 
the angular momentum quantum number n, as shown 
in Fig|5] There is a simple linear relationship E oc n 
for small \n\ (up to n = 10) when fi ^> Aq. This linear 
relationship could be explained by quasi-classical analysis 
in the next section. As ~ 0, the bound state energy 
except the zero mode merges into the continuum Ao . The 
prediction of the bound state energies by linearizing the 
vortex amplitude [2(1 [2l| is proportional to V^V (N is an 
integer quantum number), which is inconsistent with our 
calculation. 

We also solve for the bound state energies of the 
Schrodinger Hamiltonian for comparison. Fig|3] shows 
the bound state energies as a function of /i, with 
mc 2 /A = 10. Due to the time-reversal symmetry, 
= —E n _x. Unlike the Dirac case, no particle-hole 
symmetry remains. For highly positive chemical poten- 
tial jj, ^> Ao, the qualitative behavior is the same as the 
Dirac case except that the zero mode is missing. A sim- 
ple linear relationship E cx n is also recovered, as shown 
in FigJU The negative chemical potential, fi < 0, can 
be realized in superfluid Fermi gases near the BEC to 
BCS crossover [26|. For [i < 0, the continuum of the ex- 
citation becomes \/ 4- Aq. The mid-gap state merges 
into the continuum when the chemical potential is deeply 
negative, \i <S — Ao. The qualitative behavior of the mid- 
gap states remains the same for different c, which is also 
shown in FigsJS] and H] with mc 2 / 1 Aq = 1. 
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Under the external magnetic field perpendicular to 
the surface, the vortex bound states of the Dirac and 
Schrodinger Hamiltonian behave very differently. 

For the Dirac Hamiltonian, the continuum thresh- 
old of bound state energy is Ao yjl - h 2 /A 2 when 
< h < fi/y/l + A 2 ,/// 2 , and y/fi 2 + A 2 , - h when 
/t/\A + Aq//^ 2 < h < \/ n 2 + A 2 , and the energy be- 
comes fully continuum when h > y/ fi 2 + A 2 ,. 

In the limit of \n\ ^> Ao, as illustrated in Fig J7] for 
///Ao = 10, the vortex bound states of different angular 
momenta are equally spaced. As the external magnetic 
field increases, the vortex bound states merge into the 
continuum. Notice that all the states are insensitive un- 
der the external magnetic field. 

As |/t| decreases, the vortex bound states of finite an- 
gular momentum moves towards the continuum, as il- 
lustrated in FigE] for /t/A = 1. When |/t| < A , 
the zero mode exists until it meets the continuum when 
h > y/ /i 2 + Aq. The remaining bound states are located 
near the edge of continuum. The case that ///Ao = is 
shown in FigO 

For the Schrodinger Hamiltonian, it is easy to see from 
Eqs. (J5j)-(|6j) that the vortex bound state energy shift due 
to the external magnetic field is ±/i. This is the funda- 
mental difference from the Dirac Hamiltonian. 



IV. QUASI-CLASSICAL ANALYSIS 
A. Dirac Hamiltonian 

From Eq.([2]), the eigenvector of the zero energy state 
is given by 



vpp ■ a — /i — ha z 



-ia y A* 



ia y A 
-VfP- a + M • 



ho. 







(9) 

Under the condition of the quasi-classical approxima- 
tion, i.e., ^ 1, the wavelength of the quasi-particles 
is much smaller than the coherent length of the vortex, 
their trajectories are almost straight lines. Following the 
argument by Volovik [27| . the zero energy state corre- 
sponds to the trajectory that crosses the center of the 
vortex, i.e., zero impact parameter, b — 0. The super- 
conducting gap function 



A(s) = A(| S |) S gn( S ) 



(10) 



where s parameterizes the trajectory across the vortex 
center. 

For \h\ < by directly solving Eqs. ([§))- (ITU)), we have 
the zero energy state 
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where the plane wave momentum q = ^ ^ - (e x cos ( 



e y sin 9) , and 



m = cx P (-^ J 8 ds'A(\s'\) S gn(s': 



(12) 



The non-zero angular momentum vortex bound state 
near the zero energy state can be estimated by small im- 
pact parameter b, A(r = yj s 2 + 6 2 )e^ ~ A(|s|)sgn(s) + 
iA(\s\)b/s. Hence we have the perturbation 



H' 



iA(\s\)b ( ia 



itjy 



(13) 



The bound state energies up to first order correction be- 
comes 

E = (H ) = ag n(/i)6 J -~ ' ; / (14) 

and the angular momentum L z = —pb = — b\(j,\/vF which 
is nh by the semi-classical quantization condition, 



E = —nhujo 



with 



w 



(15) 



(16) 



Notice that luq has no h dependence. 

On the other hand, when \h\ > |/t|, the zero energy 
state is a decaying solution, 



s /(g)ci 



2-v/mG"- ft) 



/ h-j* \ 
—i\]h 2 — /t 2 
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(17) 



where k 



The correction (H 1 ) — 0. It is 



consistent with all finite angular momentum bound states 
merging to zero bound state at h = fx. 



B. Schrodinger Hamiltonian 

For the case of the Schrodinger Hamiltonian with 
superconducting gap, the eigenproblem for the quasi- 
particles is 
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2m 



A 

A* 



- h 



= E 



(18) 



In the quasi-classical approach, after the transforma- 
tion (u^,i>i) — ¥ e iq ' r (u^, vi), the eigenproblem becomes 
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A* 
fi 2 a. 
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- h 



iA(\s\)b 



= E 



(19) 
(20) 



where q = ^ 2 ™n ^ CQg 



+ c y sin f 



, and z% = hq/m, for 
A* > 0. 

We solve the left hand side, and treat the right hand 
side to the first order perturbation 22j . The bound state 
energy is 



E = — (n — —)foojo — h 



(21) 



with wo defined in Eq. (|16p . n is an integer. The half 
integer i is due to the contribution of the term — \t "\ Sn 
which is absent in the Dirac Hamiltonian. Because of the 
extra half integer, there is no zero mode in the vortex 
bound state of the Schrodinger Hamiltonian. 

The bound state energy for (u^,v^) can be obtained 
from the particle-hole transformation, which gives 



E = (n- -)hui + h (22) 



V. CONCLUSION 

In conclusion, we solve numerically and perform the 
quasi-classical analysis for the vortex bound states of 
both the Dirac and Schrodinger Hamiltonians in the pres- 
ence of superconducting gap. It was found that the 
bound state energies follows the linear relationship with 
the angular momentum for both cases when the chem- 
ical potential is highly positive, fj, Aq. If the chem- 
ical potential is highly negative, fi <C — Ao, the linear 
relationship holds for the Dirac Hamiltonian due to the 
particle-hole symmetry; whereas the bound states for the 
Schrodinger Hamiltonian merge into the continuum. Af- 
ter applying the external magnetic field perpendicular 
to the interface, the bound states energies for the Dirac 
Hamiltonian are insensitive; whereas there is a Zeeman 
shift for the Schrodinger Hamiltonian. The bound states 
should be determined by the local density of states on 
the surface from the scanning tunneling microscopy ex- 
periments. The bound state energy differences can be 
determined by the absorption of circular polarized light. 
The qualitative differences between the Schrodinger and 
Dirac Hamiltonian could serve as an indirect evidence of 
the existence of Major ana fermions. 
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FIG. 2: Vortex bound state energies of the Dirac Hamiltonian 
as a function of angular momentum quantum numbers n in 
the absence of external magnetic field, h = 0, for n/Ao = 
10,1,0,-1,-10. 



FIG. 1: Vortex bound state energies of the Dirac Hamiltonian 
as a function of chemical potential /i in the absence of external 
magnetic field, h = 0, for different angular momentum label 
by quantum number n. E— n = — E n . 
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FIG. 3: Vortex bound state energies of the Schrodingcr 
Hamiltonian as a function of chemical potential n in the ab- 
sence of external magnetic field, h = 0, for different angu- 
lar momentum label by quantum number n. mc 2 /Ao = 10. 
E~ n — —E n -i. The dashed lines represent the continuum 
threshold. 
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FIG. 4: Vortex bound state energies of the Schrodinger 
Hamiltonian as a function of angular momentum quantum 
numbers n in the absence of external magnetic field, h = 0, 
for n/Ao = 5, 1,0, -0.5, -1 when mc 2 /A = 10. 
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FIG. 5: Vortex bound state energies of the Schrodinger 
Hamiltonian as a function of chemical potential n in the ab- 
sence of external magnetic field, h = 0, for different angu- 
lar momentum label by quantum number n. mc 2 /Ao = 1. 
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-i. The dashed lines represent the continuum 
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FIG. 6: Vortex bound state energies of the Schrodinger 
Hamiltonian as a function of angular momentum quantum 
numbers n in the absence of external magnetic field, h = 0, 
for ^/A = 5, 1,0, -0.5, -1 when mc 2 /A = 1. 
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FIG. 7: Vortex bound state energies of the Dirac Hamiltonian 
as a function of external magnetic field h at n/Ao = 10, 
for different n. The dashed lines represent the continuum 
threshold. 
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FIG. 8: Vortex bound state energies of the Dirac Hamilto- 
nian as a function of external magnetic field h at /j/Ao = 1, 
for different n. The dashed lines represent the continuum 
threshold. 



FIG. 9: Vortex bound state energies of the Dirac Hamilto- 
nian as a function of external magnetic field h at fi/Ao — 0, 
for different n. The dashed lines represent the continuum 
threshold. 



